m 



o 



CKM matrix element Vcb in a QCD Potential model 

^'^D K Choudhury and ^Krishna Kingkar Pathak 

^Deptt. of Physics, Gauhati University, Guwahati-781014, India 

^ Physics Academy of North-East, Guwahati-781014, India 

^Deptt. of Physics, Arya Vidyapeeth College, Guwahati- 781016, India 



^ ; e-niail:kkingkar@gniail.coni 

u • 

<'■ 

Abstract 

vn ■ 

CN I We study the slope p^ and curvature C of Isgur Wise function for the heavy-hght mesons 

,_i' in general with particular emphasis on B meson in a QCD Potential model. The IW function 

i^ . is then used to compute the partial decay width and branching ratio for the semileptonic decay 

I I of (i?° -^ D,D*lu). The computed value of the CKM element \Vcb\ = 0.040 is found to be in 

Mh' agreement with the available data. 
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^ : 1 Introduction 

00 ; 

\C) ' In a recent communication [1] , we have reported the results of pseudoscalar decay constants of heavy- 

■^ \ light mesons in a QCD potential model with linear part of the potential as perturbation. The technique 
used was the quantum mechanical perturbation theory with plausible relativistic correction. 
In this paper, we extend the QCD Potential Model with its input parameters towards the semelep- 
^ ' tonic decays to study the Isgur- Wise function. Exclusive semileptonic decays of hadrons containing 

k>( ■ a bottom quark provide a path to measure the Cabibbo-Kobayashi-Maskawa(CKM) matrix elements 
;_( ! V^fe,an important parameter to test the Standard Model. It is well-known in the literature that in 

- - -' case of heavy to heavy transitions like 6 — )■ c decays, all heavy quark bilinear current matrix elements 
are described in terms of only one form factor, which is called the Isgur- Wise (IW) function in leading 
order. The IW function, particularly its slope (^'(l))at the zero recoil point is important since it 
allows a model independent way to determine the CKM element \Vcb\ for the semileptonic decay of 
[B^ — )• D*lu) and {B^ — )■ Dlv). This method of obtaining the CKM element within the framework of 
Heavy Quark Effective Theory (HQET) was proposed by Neubert|2]. He observed that the zero recoil 
point is suitable for the extraction of CKM element \Vcb\- The method basically relies on the existence 
of one universal form factor(Isgur-Wise function) and the fact that the form factor is unaffected from 
l/rnq {Q = b, c) corrections at zero recoil. 

Here, we compute the slope p^ and curvature C of the Isgur- Wise funcion for B meson and then 
obtain the CKM element Vcb- In section 2, we discuss the formalism and summarise the results and 
conclusion in section 3. 



2 Formalism 

2.1 Wave function in the model 

The wavefunction obtained by using the Dalgarno method of perturbation with hnear part of the 
Cornell Potential as perturbation is being discussed in ref.[Tl El S]. For completeness we write the 
wavefunction here again. 
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The QCD potential is taken as 
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V^(r) = Os + br + c (7) 
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Here Aq is the undetermined factor appearing in the series solution of the Schrodinger equation. However 

with Ao=0 the effect of c vanishes entirely from the solution. 

The corresponding normalised wavefunction in momentum space is[T]: 
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2.2 The Isgur-Wise function and Semileptonic decay 



In case of weak decay like the semileptonic B — > Dlu decay, the amplitude is a product of two matrix 
elements, the hadronic weak currentif^ and the leptonic weak currentL^ with its vector V^ and axial 
vector A^ i.e 

A = ^V,,{rih)Mk2)\L,\0) X {DiP')\H^\B{p)) (9) 



where H^ = C7''(l - 75)6 and L^ = l'y''{l - '^^)i'i. 

Since both the B and D mesons are pseudoscalar (J^ = 0~), so from the consideration of lorentz 
covariance and parity one has {D{P')\A^\B{p)) = and there remains the vector part 

{D{p')\V^\B{p)) = U{q'){p + p')^ + f4q'){p-p')^. (10) 



where f+{q^) and/_(g^) are two form factors. 

Similarly for B — )■ D*!^, we get another four independent form factors 

{D*{p\e)\urb\B{p)) = 2^6^-^ /-^"^^ V{q') (11) 
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Unlike the case of electromagnetic interaction here the normalisation of weak form factors are in 
general unknown. 

However, in the limit of infinitely heavy quark masses Mb — >■ oo, a new heavy flavour symmetry 
appears in the effective Lagrangian of the standard model which provides the model independent 
normalisation of the weak form factors in the framework of HQET. 

In the heavy-quark limit the masses of the heavy quarks and consequently, the masses of the heavy 
hadrons are taken to be infinite. This leads to an additional symmetry which is known as the heavy 
flavour symmetry. With the hadron masses their momenta also go to infinity but the hadron four 
velocities remain finite in this symmetry. One is then interested in the dependence of form factors on 
the (finite) velocity product v.v' . Moreover, the heavy quark symmetry is an approximate symmetry 
and correction arises since the quark masses are not infinite. This correction may be studied system- 
atically in the framework of HQET. The leading symmetry-breaking corrections are from the terms 
of the order of -^ .Although the relative corrections can be calculated using perturbative QCD, the 
VTiQ* corrections produce new uncalculable functions. For B — )■ D*lv decay there are found to be 
four such uncalculable functions. Hence the predictive power of the theory is reduced. Isgur, Wise, 
Georgi and others[2l [5l [6] showed that in case of weak semileptonic decays of -B — )■ Dlv or _B — )■ D*/i/, 
all the form factors are expressible in terms of a single universal function of velocity transfer, which 
is normalized to unity at zero recoil. This universal function is known as the Isgur- Wise function, 
which measures the overlap of the wave functions of the light degrees of freedom in the initial and 
final mesons moving with velocities v and v' respectively. 



2.3 Isgur- Wise function in the model 

In the heavy quark limit close to y = 1, the Isgur- Wise function is written as : 

eK-<) = e(y) 

= l-p^(y-l) + C{y-lf + ... (13) 

where 

y = Vf^.v'^ (14) 



and f^ and v'^ being the four velocity of the heavy meson before and after the decay. The quantity p 
is the slope of I-W function aX y = 1 and known as charge radius : 
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The second order derivative is the curvature of the I-W function known as convexity parameter : 
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For the heavy-hght flavor mesons the I-W function can also be written as [U [H] : 
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where 



/ = 2/i2 {y - 1) 
Upon integration of eqnJTTl one obtains|ll [12] 
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The normalized wave function naturally follows the zero recoil condition of the IW function ^(1) = 1 
in the model. 



2.4 The strong coupling constant ag in the Model 

In the perturbation procedure, the convergence point of view demands from the eqnj8]that 
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This condition provides a lower limit of a^ > 0.38 in the model to incorporate a small value of 
Po{Po < A^)- Again the reality condition of e provides an upper limit of a^ < |. Hence the condition 



of equation. fl20|) allows a range of 0.38 < as < 0.75 in the model to treat the linear part of the 

Potential as perturbation with its model parameters. 

Instead of using a free strong coupling constant a^ within this range, we therefore consider the strong 



running coupling constant appeared in the potential V{r), is related to the quark mass parameter 
asP l [3 [8] 

^s (/i^) = 7 ..^^'' /,..,...x (21) 
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where, nj is the number of flavours, fi is the renormalisation scale related to the constituent quark 
masses as u = 2 "^'"^^ and Mb is the background mass. The background mass is not an arbitrary 
parameter and can be calculated in the frame work of lattice QCD. The apperance of mass Mb in 
eqnJH] is similar to the case of QED where a has the mass of e^e~ pair under logarithm [7] . Here we 
relate the background mass to the conflnement term of the potential as Mb = 2.24 x 6^/^ = 0.95GeV 

mm- 

In studying leptonic decay [1] , the value of Aqcd = 200MeV was flxed in the eqnj2T] to obtain the 
value of running background coupling a^. In the case of leptonic decay of charged mesons, the quark 
and antiquark annihilate to produce a virtual W^ boson such that g^ = M^ and hence we get only 
one form factor, the decay constant fp which absorb all the strong interaction effects. In the case 
of semileptonic decay, however the case is different since q^ is different for event to event and hence 
more than one form factor appears. This decrease of q^ in semileptonic decay leads us to consider a 
larger value of Aqcd in compared to that of the leptonic decay which effectively increases the strong 
coupling constant ^^(/i). The physically plaussible range of effective Aqcd should in principle be 
deduced from the allowed bands of the slope and curvature of the I-W function. Considering the 
theoretical bounds on slope 3/4 < p^ < 1.51 |9l |T0] and curvature C > -^[lOj, we draw a curve for 
the variation of I-W function for B meson and is shown in flg.l. The allowed range of p^ provides a 
range of Aqcd in the model as 382MeV < Aqcd < ^SOMeV. 

In ref.[7j, it was analysed in studying the freezing of QCD coupling effects that for running back- 
ground coupling in V scheme, one can choose Ay{nf = 3) = 410MeP^(eqn.30 of ref. [7]) and this 
value of Ay does not contradict with those which are commonly used in MS renormalisation scheme 
and give rise to as{Mz) = 0.118 ± 0.001. Hence we justify ourself to choose Aqcd = OAlOGeV in 
this work. The input parameters used in the numerical calculation are the same as is used in our 
previous works[lJ Uf = 3, ma = 0.336GeV, rric = 1.55GeV, nib = 4:.97GeV with b = 0.183 GeV^ and 
cAo=lGe1/2/3 with c = -OAGeV 



3 Results 



3.1 Determination of CKM element Vcb 

Using eqnll3] and eqnJT9l we compute the slope and curvature of the I-W function for B meson. It 
is also added that the result with Aqcd = 410 MeV is found to be improved than our previous work 
with V Scheme [T2]. 



Table 1: Slope and curvature of I-W function for B and D mesons. 
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Figure 1: Variation of I-W function with Y for different scales of A. 



From the tableJH we see that with Aqcd = 200 MeV, the results overshoots all other theoretical 
upper bounds 3/4 < p^ < 1.510 [10] and hence for determination of Vcb we consider the value of p^ 
and Cb with Aqcd = 410MeV^. 

The differential semileptonic decay rate B ^^ Dlv for the massless leptons is given by [21 [20] 



dV G\ 



dy 487r3 



iKbP Ml {y' - If /'{Mb + Mof e {v) , 
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where y lies in the range 1 <y < ^^m m "'' • 

The differential semileptonic decay rate i? — t- D* / z/ is defined by 
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y + 1 {Mb-Md')\ 
where q'iy) = M^ + M|,* — 2yM£,*MB- By integrating the expressions for the differential decay rates 
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Table 2: Comparison of slope and curvature of B mesons with other works. 



work 


'2 

Pb 


Cb 


Present 


0.993 


1.114 


Faustov etal; 13 


1.04 


1.36 


Lattice QCDfH] 


o.sstlltlt 




ALEPH 15 


0.92 ±0.98 ±0.36 




Belle [16] 


1.12 ±0.22 ±0.14 




Le Youanc et al (10 


>0.75 


>0.47 


QCD Sum Rule [T7| 


0.65 


0.47 


Relativistic Three Quark Model [18] 


1.35 


1.75 


Neubert [19] 


0.82±0.09 





1± 



Ay 



q^{y) 
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y + l {Mb-Md*Y 
we get the predictions for the total decay rates in our model as 

T{B -^ Dlu) = 6.82\Vcb\^ ps'^ 
T{B -^ D*liy) = 28A0\Vcb\'^ ps'K 



(25) 



(26) 

(27) 
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Taking mean values of lifetimes from PDG2012: r^o = 1.519 x 10~ s and tb+ = 1.641 x lO^^"' s, 
and also using the experimental masses for the eq.24 and eq.25, we compute the semileptonic branching 
ratio by using the relation BR = P x r and find 



BR{B° -^ D+ru) = 


= 10.401 Va ^ 


BR{B+ -^ D'^l+u) = 


= ii.ioiKfer, 


BR{B^ -^ D*+l-v) = 


= 43.31|K6r, 


BR{B+ -^ D*^l^u) = 


= 46.78|K6|l 



(28) 

The comparison of these theoretical results with the experimental branching ratios from the lattest 
PDG [2T] gives us following values of the CKM matrix element \Vcb\'- 



BR{B^ -^ D+ruy'P = 0.0218 

BrIb+ -^ D°/+z/)^"P = 0.0229 

fii?(fi° -^ D*+l-iyy''P = 0.0509 

BrIb+ -^ D*°l+iyy''P = 0.0558 



\Vcb\ 
\Vcb\ 
\Vcb\ 
\Vcb\ 



0.046, 
0.045, 
0.034, 
0.035. 



(29) 



Thus the averaged \Vcb\ over all presented measurements of semileptonic decays B — > Deu and B — )■ 
D*ez/ is equal to 

\Vcb\ = 0.040 (30) 



and is in good agreement with the experimental result |21] . 

\Vcb\ = 0.0396 ± 0.0009 (exclusive). 



4 Conclusion 

In this work, we have studied the renormahsation scale dependance of Isgur-wise function by using 
a wavefunction with hnear part of the Cornell Potential as perturbation. Considering the exclusive 
semileptonic decay of 5 — >■ Deu and B — )• D*ev, we then compute the CKM element | V^;, | in this 
approach. The result of the CKM element is found to be within the error limits of other results. The 
following features are observed in this work: 

• The renormalization scale of the model was set to be Aqcd = 410 MeV as is used in ref. [7]. 

• The slope and curvature of the Isgur wise function are found to lie within the range of limits 
found in the literature. 

• The computed value of the CKM element from the exlusive semileptonic decay of B meson is 
obtained as \Vcb\ = 0.040. The result is to be well consistent with the lattice result \Vcb\ = 
0.0409 ± 0.0015 ±0.0007[ 



• It becomes an worthy comment to note that the result with Aqcd = 200 (as is done in our 
previous work) is found to be very poor to study the branching ratio and \Vcb\- 

From a phenomenological point of view, the present theoretical framework is a considerably useful 
tool to investigate various physical quantities for the heavy-light quark systems. Further work on the 
renormahsation schemes and scale dependance of strong coupling as is under study. 
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